Energy shifts of Rydberg atoms due to patch fields near metal surfaces 
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The statistical properties of patch electric fields due to a polycrystalline metal surface are calcu- 
lated. The fluctuations in the electric field scale like 1/z 2 , when z S> w, where z is the distance to 
the surface, and w is the characteristic length scale of the surface patches. For typical thermally 
evaporated gold surfaces these field fluctuations are comparable to the image field of an elementary 
charge, and scale in the same way with distance to the surface. Expressions for calculating the 
statistics of the inhomogeneous broadening of Rydberg-atom energies due to patch electric fields are 
presented. Spatial variations in the patch fields over the Rydberg orbit are found to be insignificant. 

PACS numbers: 34.35. +a, 32.80.Ee 
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I. INTRODUCTION 

Excitation to a Rydberg state enhances an atom's in- 
teraction with a metal surface. At large atom-surface 
distances, this results in energy level shifts that can be 
calculated using first-order perturbation theory At 
smaller distances, the influence of the surface is more 
drastic — the Rydberg atom can be "field-ionized" by the 
surface @, 0). These phenomena may be visualized as 
arising from the interaction of the Rydberg atom with 
the electric fields due to its electrostatic "image." Com- 
pared to an atom in the ground state, a Rydberg atom 
has an enhanced susceptibility to these fields. This is be- 
cause the Rydberg electron experiences a greatly reduced 
electric field from the ion core due to their larger average 
separation. 

Polycrystalline metal surfaces generate inhomogeneous 
"patch" electric fields outside of their surfaces [4J . These 
fields may also influence Rydberg atoms, potentially 
causing both level shifts and ionization and competing 
with the more intrinsic image charge effects. In general, 
patch fields arise from the individual grains of a poly- 
crystalline surface exposing different faces of the bulk 
crystal. Each face has a different work function due to 
differing surface dipole layers @. For example, Singh- 
Miller and Marzari [6( have recently calculated the work 
functions of the (111), (100), and (110) surfaces of gold 
and found 5.15, 5.10, and 5.04 eV, respectively. These 
differing work functions correspond to potential differ- 
ences just outside the surface beyond the dipole layer. 
Consequently, charge density must be redistributed on 
the surface to satisfy the electrostatic boundary condi- 
tions, producing macroscopic electric fields @. While 
patch fields were first discussed extensively in the con- 
text of thermionic emission [4( , they may be present near 
polycrystalline metal structures of any type, including 
electrodes and electrostatic shields. 

Recent advances in the trapping of cold atoms near 
surfaces have opened up the possibility of precision stud- 
ies of Rydberg-atom-surface interactions as a function 
of atom-surface distance. For example, Tauschinsky et 
al. 0| have recently observed electromagnetically induced 



transparency due to Rydberg excitation of atoms at 
10 fim to 200 ^m away from a gold surface. Lesanovsky et 
al. § have calculated some interesting properties of Ry- 
dberg atoms exposed to inhomogeneous magnetic fields 
due to magnetic microtraps, and Crosse et al. Q have 
recently calculated level shifts and transition rates of ru- 
bidium atoms near a copper surface at room temperature. 
In addition, there are several quantum information pro- 
posals that will involve Rydberg atoms in proximity to 
metal surfaces [Toj - TT^ I . Consequently, it is desirable to be 
able to estimate the influence of patch fields on Rydberg 
atoms. In this paper, we examine relevant models of the 
surface, report on the statistics of the patch fields, and 
determine the influence of these fields on Rydberg-atom 
energies. We assume that the atom-surface distance is 
large compared to other relevant length scales and that 
the atomic energy level shifts can be treated using per- 
turbation theory. 

The model we adopt for the patch fields is similar 
to one used by Rzchowski and Henderson Their 
work was motivated by the Witteborn-Fairbank experi- 
ment which was intended to compare the force of 
gravity on electrons and positrons. Due to the rela- 
tively weak gravitational force, electrostatic shielding was 
necessary — the charged particles traveled down the axis 
of a hollow copper cylinder used for shielding. It was 
important to understand the variations in electrostatic 
potential along the axis of this tube due to patch fields, 
and Rzchowski and Henderson obtained results relevant 
to this geometry. In the present work, we concentrate 
on a planar surface and the statistical properties of the 
electric field and its spatial derivatives. 



II. RYDBERG-ATOM ENERGY SHIFTS IN 
EXTERNAL FIELDS 

We first calculate the energy-level shifts of a sin- 
gle atom in response to the local electrostatic potential 
V(x,y,z) created by the patches. This will allow us to 
calculate the statistics of the energy shifts once the statis- 
tics of the patch fields are known. 
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We consider an addition to the atom's Hamiltonian Hq 
of the form [H 



i 



(1) 



where pi and Qij are the dipole and quadrupole moment 
operators, respectively, and Di is the operator represent- 
ing the derivative with respect to the ith argument. The 
quadrupole and higher-order moments will allow us to 
consider the influence of an electric field varying over the 
extent of the atom — of particular interest because Ryd- 
berg atoms are much larger than ground-state atoms. 

We will consider the shift given in Eq. (fTJ) using first- 
order perturbation theory. In the absence of preexisting 
external fields, there would normally be no contribution 
from the dipole term to the first-order shift. However, 
we will assume that a dc electric field aligned with the 
surface normal z has been applied. This may be done to 
enhance sensitivity to patch fields, break degeneracies, or 
for technical reasons (see, for example, Ref. [12J). The ef- 
fect of this field is incorporated into H . Our basis states 
(eigenstates of H ) will be considered to have a cylin- 
drically symmetric charge distribution about the surface 
normal. This symmetry constricts the moments, so that 
f^x l^y — 0; Qxx — Qyy — "^Qzzi etc. 

Equation (fTJ) involves the evaluation of arbitrary order 
mixed derivatives of the potential in all three spatial di- 
mensions. However, the introduction of cylindrical sym- 
metry allows a considerable relaxation in this require- 
ment if spherical (instead of cartesian) multipole mo- 
ments are used. Therefore, the external potential due 
to the patches is expanded in the form 



V 



(2) 



about the location of the atom, where dg^ m are the ex- 
pansion coefficients, r is the distance away from the 
center of the expansion, 9 and <f> the normal spherical 
coordinates with the polar axis aligned with the sur- 
face normal, and Cg m are rescaled spherical harmonics; 

Ct, m = y/An/(2e+l)Y t , m . 

To obtain an expression analogous to Eq. (TTJ) , this new 
expansion of the potential is substituted into the volume 
integral for the electrostatic energy due to a charge dis- 
tribution p in an external potential, E — J dr pV, with 
dr as the differential volume element. We obtain 



E = 2_, dt,o dr p r e P e (cos ( 

B J 



(3) 



where terms involving m ^ do not appear, due to the 
cylindrical symmetry of the charge distribution, and Pg 
are the Legendre polynomials (= Cf.fi). The values of 
the dg.fi coefficients in the expansion can be readily de- 
termined from Eq. ^ by evaluating the derivatives of 



the potential with respect to the distance to the sur- 
face z evaluated at the origin of the expansion, giving 
df.fi = (1 / £1) D l z V (x , y , z) , where D £ z means take the £th 
derivative with respect to the z coordinate. The first 
order energy shift can be written in a form that only 
depends on the gradients of the field in the z direction 
evaluated at the location of the atom: 



E 1 = YM e [DiV(x > y,z)] ) 



(4) 



where Mg has been introduced to simplify notation. To 
evaluate these matrix elements, we assume that only the 
charge distribution due to the Rydberg electron needs to 
be accounted for, so that Mg = (l/£!)g e (^>o|^C'e,o|'0o}> 
where q e is the electron charge, and |^o) are the energy 
eigenstates of the zeroth order Hamiltonian Hq. The val- 
ues of Mg are proportional to the normal spherical mul- 
tipole moments (see, for example, Ref. (16|). 



III. STATISTICS OF THE PATCH FIELDS 

As shown in the previous section, the energy of any 
particular atom depends on the field at its location. Con- 
sider an ensemble of atoms placed a certain distance z 
away from the surface. In general, the patch fields are 
statistical in nature, so that spatial inhomogeneities in 
the field will cause an inhomogeneous broadening in the 
ensemble. We can characterize this by the variance in 
the energy of a given state, calculated using Eq. (gj), as- 
suming that the average shift is zero: 



((AE) 2 )c =Y^M t M t 



[D e z V(x,y,z)\ [D e z V(x,y,z) 



(5) 



where (--^c is used to specify an ensemble (classical) 
expectation value. 

Therefore, to calculate the variance of atomic en- 
ergy levels ((AE) 2 )c due to the statistical fluctuations 
in the field above the surface, we will develop expres- 
sions for ([D z V(x,y,z)][D z V(x 7 y,z)})c- For example, 
the most important statistical fluctuation for Rydberg 
energy level shifts is of the electric field in the z di- 
rection, which can be characterized by its root mean 
square (rms) value: \{F Z ) — (i 7 ^) 2 ] 1 / 2 , which is given by 
{([DlV(x,y : z)][DlV(x : y,z)}) c }^. 

To calculate these statistical averages, we start by con- 
sidering the solution of Laplace's equation W 2 V(x, y, z) — 
above a plane surface when the potential on the surface 
is specified. One particular solution of Laplace's equation 

Vq e ih xX+ik y y e -kz ^ w h ere £ = 



is: 



V(x,y,z) 



and k x , k y and Vq are constants. Consider the following 
superposition of similar solutions (all integrations are as- 
sumed to run from negative to positive infinity, unless 
otherwise specified): 



V(x, y,z)= / dk x dk y V(k x , k y y k * x+lk y y t 



-kz 



(6) 
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We may use this expression to determine the potential 
over any surface in the plane z = with a defined po- 
tential V s (x, y) by using the inverse Fourier transform to 
determine V(k x , k y ): 



V(k x , ky) 



1 



dxdy V s (x, y)e 



(2<p 

Putting Eqs. ((6|) and (O together gives 
1 



(7) 



2 / ""-a:""*? 



(27T) 

x ( dx'dy'V s (x',y')e- tk * x '- ik y y '. (8) 



V(x,y,z) 



Consider the covariance between derivatives of the field 
evaluated at two points a and b in space, determined 
using the proceeding equation: 



D p V(x a ,y a ,z a ) D q V(x b ,y b ,z b ) 



dk x ^ a dky }a dk X}b dky >b dx' a dy' a dx' b dy b 



cxp\ih x a x a -h iky_ a y a k a z a ik x a x a ^ky^ a y a 

x(a p , a y(al b f x (V s (x' a ,y' a )V s (x' b ,y' b )) c , (9) 

where a\, a = ik x ,a, ot2, a = ik y , a , Q3,a = —k z , a , ai, b = 
ik x ,b, &2,b = ik Vl b and a^j = —k z , b . As Eq. J5| shows, 
we only need p = q = 3 (derivatives in the z direction), 
but it is not difficult to deal with this slightly more gen- 
eral form, which will also allow us to calculate additional 
quantities, possibly of use to others, such as the total rms 
electric field. 

We now make an assumption about the statistical na- 
ture of the field: the correlation function C for the surface 
potential only depends on the separation between the two 
points a' and b' (i.e., it is a "stationary" process): 

(V s (x' a , y' a )V s (x' b , y' b )) c = C(x' b - x' a ,y' b - y' a ), (10) 

and rewrite Eq. (0) using Ax' = x' b — x' a , and Ay' = 

v'b - y'a- 



D e p V(x a ,y a ,z a ) D e q V(x b ,y b ,z b ) 



dk Xi adky,adk Xtb dk y , b dx' a dy' a d(Ax')d(Ay') 



(27r) 4 J """ x > a y> a x > b y>' 

e~xjp\ik x a x a -\- ik y a y a k a z a ik xa x a iky^ a y a 

ik x fiX b iky^ b y b k b z b 
+ik x , b (Ax' + x' a ) + iky, b (Ay' + y' a )} 
x(a p , a ) £ « b ) £ ' xC(Ax',Ay'). (11) 

Use of the familiar relationship: 

/ dx' exp[ix'(k a — k b )] — 2ir8{k a — k b ), where #(•••) is the 



Dirac 5 function, allows simplification to 



D p V(x a ,y a ,z a ) D q V(x b , y b , z b ) 



dk x dky(a P Y(a*Y 



(2tt) 2 J —y^p> 

x exp[ik x (x a - x b ) + ik y (y a - y b ) - k(z a + z b )] 



x / d(Ax')d(Ay')exp[ik x (Ax') + iky(Ay')} 
xC(Ax',Ay'). (12) 

Assuming the surface has no preferred direction, 
C(Ax', Ay') is only a function of Ar' = V 'Ax' 2 + Ay' 2 , 
and the evaluation of the last integral in Eq. [12] is equiv- 
alent to taking the two-dimensional (2D) Fourier trans- 
form of a radially symmetric function (see, for example, 
Ref. [13): 



W(k) 



1 



d(Ax')d(Ay') exp[ifc x (Aa;') + ik y (Ay')} 



2ir 

xC(Ax',Ay') 

d(Ar')Ar'J {kAr')C(Ar'), 



(13) 



where Jo(- • • ) is the zeroth-order Bessel function. Equa- 
tion ([T2|) may then be written as 



D e p V(x,y,z) D e q V(x + Ax,y + Ay,z + Az) 



= —j dk x dk y (a P Y(a* q Y W(k) 

x exp[ik x Ax + iky Ay — 2kz — kAz]. 



(14) 



A generalisation of this result to mixed derivatives is 
straightforward, but the notation cumbersome. To eval- 
uate Eq. ([5]), we need a slightly less general expression: 



(^D l z V{x,y,z)\ [D z 'v(x,y,z)]) 



(-1) 



e+e' 



dk W(k)k 



1+1 



exp[— 2k z] 



(15) 



It is helpful to rewrite this in a dimensionless form. 
A natural length scale for the surface is w = l/Vd, 
where d is the mean areal density of the surface patches. 
We assume that the covariance of the surface poten- 
tial depends on w in such a way that it can be writ- 
ten in terms of a scaled covariance function C as 
C(Ar') = $ 2 ms C*(Ar'/ W ), where $ rras ~ [{V, (x, y) 2 ) - 
(V s (x, y)) 2 ] 1 ' 2 is the rms variation of the surface poten- 
tial from the mean. We now introduce 



W(u) 



d( Ar' /w) ( Ar' /w) J (u Ar'/w)C( Ar' /w) , 



(16) 
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which allows us to rewrite Eq. (|15l) as 

([DfV(z,y,*)] [l>ji'V(a;,y,z)]) 



= (-lY +£ ' / du W(u) x u 1+e+t x exp[-2 M (z/w)}. 
Jo 

(17) 

In general, for the W(u) that we are interested in (see 
below), these integrals do not have closed forms. How- 
ever, they may be approximated for large z/w using an 
asymptotic technique. Part of the integrand, W(u) x 
yi+t+l ^ ma y De wr itten as a Taylor series in u about 
u = 0. Once multiplied with the rest of the integrand 
(exp[— 2u (z/w)}), the terms in the resulting series can be 
individually integrated in closed form (see, for example, 
Ref. [18|). Introducing G(l + £', z/w) as a shorthand for 
the left-hand side of Eq. (fT7)) . we obtain: 



G(L,z/w) = 

(-d l E 



i=0,2,4... 



«•<■»(„, , (») 



L+2+i 



(18) 



where 1/ = t + £' and W'W(O) is the ith derivative of 
W(u) evaluated at u = 0. Note that from its definition 
[Eq. (fTB"|) ]. the odd derivatives of W"(u) vanish at u = 0. 
For use later in this paper, we write out the first few 
terms of this series for small L: 



(19a) 



G(i,^)=-^(o)(f) a -|^(o)(H) 6 + 



G(2, zM J-W {Q )^)\ 1 ^(0)^) 6 + 

G(Sl ,/»)-|^(o)(H) 5 _«^)(o,(H) r - 



(19b) 



(19c) 



(19d) 



The first terms of these series are almost certain to 
dominate when z 3> u>. From the definition of W(u) in 
Eq. [T6l it can be seen that W(0) = requires that the 
covariance function satisfies C(Ar') < over some range 
of Ar', so that the integral taken over Ar' is zero. This 
can be interpreted physically as antiferroelectric ordering 
of the surface potential; a case which seems unlikely to 
apply to polycrystalline metal surfaces. 

It is important to note that, subject to the assumptions 
above, the details of C(Ar') do not affect the (z/w) scal- 
ing of G(L, z/w) but only its magnitude. Therefore, the 
z scaling of the patch fields is independent of the form of 
C(Ar'). 



IV. MODELS FOR THE SURFACE PATCH 
POTENTIALS 

We will now calculate G(L, z/w) using several different 
models for the electrostatic potential distribution on the 
surface. We start by calculating C(Ar') for the model 
and then use this to find W(0) and thus G(L, z/w). 

A commonly used model for the surface potential co- 
variance is of the form [l3|, l20j 



C(Ar') 



^ rms° 



(20) 



where 7 is dimensionless and on the order of 1. This 
model follows from Poisson waiting statistics for grain 
boundary crossings. This, however, is an assumption, 
and a formal justification does not appear in the liter- 
ature. An advantage of this model is that W(u) has a 
closed form [using Eq. (TTT 



W(u) 



7 



[ 7 2 + u 2]3/2 ' 



(21) 



and thus the coefficients in the expansion of Eq. ([19)) 
are readily determined [W(0) = I/7 2 , W<-^(0) = -3/7 4 , 
etc.]. 

Motivated to provide a justification for Eq. ([20)) (and 
determine a specific value for 7), we performed Monte 
Carlo simulations to calculate a surface potential covari- 
ance function according to the following recipe: 1) A total 
of N patch "centers" were randomly put within a square 
with sides of length w\f~N (for a mean areal patch density 
of 1/w 2 ). 2) At the center of this square, the patch with 
the closest center was determined. 3) As we move out 
from the center of the square in a specific direction, even- 
tually another patch center becomes closer in distance 
than the initial one. The point at which this happens is 
considered to be at a grain boundary, and beyond this 
point there is zero correlation between the local poten- 
tial and the potential at the starting point in the center 
of the square. 4) By repeating this process (generating 
N new patch centers within the square, and traveling 
out from the center until a grain boundary is reached), 
we may accumulate a surface potential correlation func- 
tion. Provided N is sufficiently large, this model seems 
physically reasonable — we are assuming that grains have 
grown isotropically outwards from randomly placed cen- 
ters on a surface. Figure [T] illustrates the results of one 
of these Monte Carlo simulations. A least-squares fit to 
Eq. (E01) gives 7 « 1.9, so that W(Q) « 0.28. 

We find that, instead of Eq. (|2"0")l . the covariance is a 
better fit to the relationship: 



C(Ar') 



$ e 

rms 



-71 ( 



-72 



(22) 



with 71 rj 1.144(4) and 72 ~ 0.993(6). The covariance 
falls off faster with increasing separation in this model. 
Man et al. [Tq] compared experimentally measured co- 
variance functions with a model similar to Eq. (|20[) and 
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Ar'/w 

FIG. 1. Monte Carlo simulation of the surface potential co- 
variance function and two least-squares-fit models. In the 
lower plot a logarithmic vertical axis is used to illustrate the 
differences at large Ar'/w. 



also found that, although exponential decay was exhib- 
ited for small separations, the covariance falls off faster 
for increasing separations (see their Fig. 2). However, 
a detailed comparison with our model is not possible as 
their surface was not isotropic. 

We have tested this model by analyzing a scanning 
electron microscope (SEM) image of an evaporated gold 
structure on a silicon substrate (see Fig . [5]) . The "wa- 
tershed" segmentation algorithm [2l], [22[ was used to de- 
termine the location of the grain boundaries. To calcu- 
late the covariance function we assume that the potential 
measured at two points separated by Ar' is perfectly cor- 
related if both points are on the same grain and uncor- 
related if the points are on different grains. As Fig. [2fc) 
shows, the computed covariance is in good agreement 
with the Monte Carlo simulation [and thus also with the 
fit of Eq. ([22])]. 

Unfortunately, a closed form for W(k) does not appear 
to be possible for the model of Eq. (|22|) . Nonetheless, it 
is possible to numerically compute the W^ l \0) required 
in Eq. (TK)1) for any C using [see Eq. ([15])]: 

/>oo 

W'W(0) = [i?Wjo(0)] / da C(a)a 1+ \ (23) 

JQ 

Performing these integrations for the model of Eq. (f2"2"j) , 
we find: W(0) = W^(0) « 0.207, and W (2 \0) « 
-0.064. 
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Ar'/w 

FIG. 2. (a) SEM image of a gold surface obtained by thermal 
evaporation (this is a portion of Fig. 4 of Ref. 23J). (b) Grain 
boundaries (indicated by white lines) over a small region of 
the image as determined by watershed segmentation. The 
average area of a patch is w 2 = (44 nm) 2 . (c) Computed co- 
variance of the surface potential based on segmentation of the 
SEM image. To calculate this from the segmented image we 
assume constant surface potentials within grains, and com- 
pletely uncorrelated potentials between grains. The Monte 
Carlo simulation of Fig. fT] is also shown for comparison. 
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numerical integration 

first non-zero term 

in expansion 

first two non-zero terms 

in expansion 




z/w 



FIG. 3. Comparison between numerical calculation of 
G(2, z/w) for the model of Eq. JUJ (using Eq. (0), and the 
asymptotic series expansion of Eq. (I19c|) . with the ^^'(0) 
coefficients calculated using Eq. (I23|) . 



FLUCTUATIONS IN THE ELECTRIC FIELD 



The G(2, z/w) function determines the variance of the 
patch electric fields: (F 2 ) c « (3/8) x ($ 2 ms /w 2 ) x W(0) x 
(w/z) 4 . Using Eq. (fl4| . the variances of the x and y 
components of the electric field can be calculated. We 
find that they are each 1/2 of the result for z. Thus, 
we can summarize; with the model of Eq. (|22[) . the rms 
electric field for z ^> w is 



E X1 



0.39- 



w 



(24) 



This result is not especially sensitive to the particular 
patch model. For example, we have performed numer- 
ical simulations of the patch field over a large array of 
square patches (each to by tu) with random potentials 
at distances z 3> w and found that the numerical pref- 
actor in Eq. flU} is 0.33 instead of 0.39. The model of 
Eq. (|20[) gives a numerical prefactor of 0.46. An approx- 
imate estimate similar to Eq. (|24l) has been provided by 
Sandoghdar et al. 24 1 and used by Mozley et al. [2oT |. 

It is worth asking when the higher-order terms of 
Eq. p^|) can be neglected. In Fig. [31 we calculate 
G{2 1 z/w) by direct integration of Eq. (Ti"5)) . The results 
due to the first two terms of Eq. (|19c|) are also shown. 
The figure indicates that keeping only the first term is 
an excellent approximation for z 3> w (a similar plot for 
the model of Eq. (j2T)]) is given in Fig. 1 of Dubessy et 
al. HI). 

The rms patch field and the image field of an 
elementary-charge both scale in the same way with dis- 
tance to the surface, so it is interesting to compare 
their magnitudes. If we assume a potential fluctuation 
of $ 2 ms = (90 mV) 2 and w — 50 nm, typical of ther- 
mally evaporated gold surfaces [H, [27| . we find that 
the rms electric field due to patches is approximately 5 
times that of the elementary charge image field: \Ei\ = 
g e /[47reo(2z) 2 ]. Seeing the influence of the image field 



due to an elementary charge near such a surface would 
be difficult. 

Despite its simplicity and intrinsic importance, there 
does not appear to be any clear experimental observa- 
tions that would support the validity of Eq. (|2~41 . Initial 
experiments with Rydberg atoms using microtrap tech- 
nology have observed large dc fields due to the deposi- 
tion of Rb on the surface Q (see also Ref. [28j]), pos- 
sibly masking the influence of patch fields. Alkali ad- 
sorption has been recognized as a problem since the very 
early days of Rydberg atom-surface interaction experi- 
ments [291 ] . Some theoretical work on the influence of 
adsorbates has been done in the context of ion-surface 
collisions [3(3]. To avoid the problem of adsorbates, Dun- 
ning's group switched to using xenon Rydberg atoms in 
their surface studies 

Dunning's group has recently studied Ryd berg atom 
image field ionization using Au(lll) samples [31j. The 
surfaces consisted of multiple grains, typically 300-500 
nm in size. Possibly due to contamination, the surface 
potential was inhomogeneous, with variations of up to 
70 mV from the average over length scales of 50-250 
nm (shown in Fig. 1 of Ref. [31]). By using a scanning 
Kelvin probe to measure surface potential, they com- 
puted the statistical properties of the electric field above 
the surface and found that the field is consistent with 
their observed image field ionization. For z > 60 nm, 
they found -E rms ~ (5 x 10~ 10 Vm)/z 2 . Assuming our 
polycrystalline model is applicable, with $ rms ~ 22 mV 
(determined from our analysis of Fig. 1 of Ref. [3l|) and 
w = 100 nm, Eq. ([2~4"|) predicts slightly larger fields, with 
-E-rms ~ (9 x 10 -10 V m)/z 2 . Given the uncertainty in de- 
termining w from the figure, this is reasonably consistent 
with the result in Ref. [3l| . 

It is possible that polar or polarizable contaminants 
at grain boundaries could reduce the magnitude of the 
patch fields. Darling [32j did extensive scanning Kelvin 
probe measurements of the work function immediately 
above copper surfaces with large grain sizes, found that 
fluctuations were significantly less than one might expect, 
and attributed the reduction to oxidation of the surface 
and physisorbed molecules (e.g. water). 



VI. PATCH FIELDS AND RYDBERG ATOMS 
ESTIMATES OF ENERGY LEVEL SHIFTS 



The statistical properties of the patch fields may now 
be combined with the atomic properties to predict the 
variance in the energy levels using Eq. ([S]). Writing this 
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as a series in w/z: 

{{AE) 2 ) C « {M 1 f{Fl) c + 2M 1 M 2 (F z d z F z ) c + ... 

«(M 1 ) 2 x(3/8)^(0)($? ms / W 2 )(^) 4 

+2M X M 2 x (-3/4)W(0)(<S> 2 m Jw 3 ) (^)' 



(25) 



The first term in this expansion is due to the rms z 
field and the atom's electric dipole, and is expected to 
be dominant at large z. However, the higher order terms 
in Eq. (|25p can be enhanced relative to the lower-order 
terms by increasing n (increasing the size of the atom). 
As the classical outer turning point of the Rydberg elec- 
tron is ~ n 2 (in atomic units), the multipole moments 
of order £ scale with n like Mi ss n 21 (see, for example, 
Ref. [33|]). These higher-order multipoles sense the field 
variations over the Rydberg orbit. 

Under what conditions will variations in the patch 
fields over the extent of individual atoms contribute to the 
inhomogeneous broadening? We may estimate this by 
equating the first two terms written explicitly in Eq. (1251) . 
This tells us that the size of the Rydberg atom, n 2 , has 
to be approximately the distance of the atom to the sur- 
face before these would be comparable. Due to the in- 
teraction of the Rydberg atom with its image, this is a 
highly nonperturbative situation We conclude that 
it would be difficult to observe any effect of the varia- 
tion in patch fields over the orbits of individual Rydberg 
atoms (at least when they have dipole moments of order 
n 2 ). An additional qualitative justification is given in the 
Appendix. 

We now give a simple numerical estimate for the inho- 
mogeneous broadening of Rydberg energy levels due to 
patch fields. When z is large compared to w and the atom 
size, the first term in Eq. (|25p dominates and the rms 
broadening will be SE = M\^J (F 2 )c- For the extreme 
Stark states of hydrogen we have Mi = fx z = (3/2)n(n — 
1), which for n = 30 is 1.7GHz/(V/cm) (this is also 



reasonable for non-hydrogenic atoms, assuming a large 
enough dc field is applied) . For a typical thermally evap- 
orated gold surface [26|, [27|, we assume 3>fms = (90 mV) 2 
and w = 50 nm, giving y/ (F 2 )c ~ 0.13V/cm. We find 
that SE 200 MHz - which should be straightforward 
to observe in optical excitation. For Rb atoms, a pos- 
sible spectroscopic probe would be the last step in the 
5s — > hp — > 5g? 5 /2 — > nk excitation sequence (where the 
last transition is enabled by a dc field sufficient to mix / 
character into the reddest nk states). 

The extreme Stark states provide the largest broaden- 
ing. Broadening due to patch fields will be much lower 
than this estimate if low angular momentum states are 
excited at fields small enough so that the Stark effect is 
second order. 



VII. SUMMARY AND OUTLOOK 

Rydberg atoms with permanent electric dipole mo- 
ments have a high sensitivity to electric fields. We have 
shown that the patch fields near a typical metal surface 
can be large compared to the image field of an elemen- 
tary charge and should be expected to cause measur- 
able inhomogeneous broadening of Rydberg energy levels. 
The rms spatial variation in the field strength has a dis- 
tance dependence of 1/z 2 . Spatial variations in the fields 
over the Rydberg atom orbit do not appear to be impor- 
tant. An experiment to verify the magnitude of the rms 
field and the expected scaling with surface distance [see 
Eq. ([24|) ] would be useful in assessing the feasibility of 
coherently manipulating Rydberg atoms near polycrys- 
talline surfaces and in planning future experiments. 
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Appendix A: Potential and field covariance functions 

There is a qualitative way to understand why the in- 
homogcncitics in the patch fields over the extent of in- 
dividual Rydberg atoms would be difficult to observe. 
Equation (|f 4p can be used to determine the covariance 
between the potential and derivatives of the potential 
measured at different locations in space. Using the re- 
sult for the 2D Fourier transform of a radially symmetric 
function we obtain: 



D e z V(x,y,z) D e z V(x + Ax,y + Ay,z + Az) 



(-1) 



, <J> 2 



'j x / du W(u)u 1+t+e ' 



xJ (uAr/w) exp[— u(2z + Az)/w], 



(Al) 



where Ar = y Ax 2 + Ay 2 . 

Again, like with Eq. (fP7|) , this integral can be approx- 
imated for large z/w by writing the W(u)u 1+e+i part 
of the integrand as a Taylor series and then integrating 
the individual terms. For covariances in the potential, 
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we obtain for the first nonzero term 
(V(x, y, z)V(x + Ax,y + Ay, z + Az)) c 

1 ~ „ / m \ 2 



Az 
2z 







1 + 


( Ar ) 2 ] 


\2z+Az J 



3/2 ' 



(A2) 

For covariance in the z-component of the electric field, 
we obtain for the first nonzero term 



DlV(x,y,z) DlV(x + Ax,y + Ay,z + Az) 

_ 3 ( Ar \ 
2 ^ 22+Az J 



3IT(0)$ 



2 

rms 



/ID\ 4 1 



Az 

2z , 



Ar 
22+A2 



7/2 ' 



(A3) 



Higher order terms involve larger powers of wjz. These 
results have been written in a way to emphasize the influ- 
ence of nonzero Ar and Az as a correction factor to the 
Ar = Az = result. It is apparent that z, the distance 
to the surface, sets the length scale for spatial variations 
in the potential and fields. Thus, we can understand in a 
qualitative way the results of the main text: a Rydberg 
atom should have a size comparable to its distance from 
the surface for spatial variations to be significant. 



If an atom is moving near a surface, spatial variations 
in the fields manifest themselves as time-dependent vari- 
ations experienced in the atom's frame. In this case, we 
note that the calculations of this section could be adapted 
to determine the power spectral densities of these fluctu- 
ations (using the Wiener-Khinchin theorem). 
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